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Higgs self-coupling measurements at the LHC 
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Both the ATLAS and CMS collaborations have reported a Standard Model Higgs-like excess at 
around mt = 125 GeV. If an SM-like Higgs particle is discovered in this particular mass range, an 
important additional test of the SM electroweak symmetry breaking sector is the measurement of 
the Higgs self-interactions. We investigate the prospects of measuring the Higgs self-coupling for 
nih — 125 GeV in the dominant SM decay channels in boosted and unboosted kinematical regimes. 
We further enhance sensitivity by considering dihiggs systems recoiling against a hard jet. This 
configuration exhibits a large sensitivity to the Higgs self-coupling which can be accessed in subjet- 
based analyses. Combining our analyses allows constraints to be set on the Higgs self-coupling at 
the LHC. 



I. INTRODUCTION 



The Standard Model (SM) Higgs Q has recently been 
excluded at the 95% confidence level from 129 (127.5) 
GeV to 539 (600) GeV by measurements performed by 
ATLAS (CMS) [2|. In addition, the Higgs mass bound 
from LEP2 was raised from 114.4 GeV [J to 117.5 GeV 
by ATLAS. Both ATLAS and CMS have also observed 
tantalizing hints for a SM-like Higgs at a mass ~ 
125 GeV with local significances of 2.5er and 2.8a, re- 
spectively. In the same mass region, the D0 and CDF 
collaborations observe an excess with a local significance 
of 2.2cr for the combination of their data sets [J]. 

Breaking down these results into the individual search 
channels has triggered some effort to pin down the prop- 
erties of the observed excess in the SM and beyond Q . 
These analyses are the first steps of a spectroscopy pro- 
gram which targets the properties of a newly discovered 
particle if the excess at 125 GeV becomes statistically sig- 
nificant. Strategies to determine spin- and CV quantum 
numbers, and the couplings to fermions and gauge bosons 
of a 125 GeV resonance with SM-like cross sections have 
been discussed in the literature @, 0| • A determination 
of the Higgs self-interaction, however, which is crucial for 
a measurement of the symmetry breaking sector remains 
challenging in the context of the SM (this can change in 
BSM scenarios Even for scenarios where the Higgs is 
close to the h — > W + W~ threshold, statistics at the LHC 
in pp —> hh + X is extremely limited 0-[H| , so that for- 
mulating constraints on the Higgs self-coupling requires 
end-of-LHC-lifetime statistics if possible at all* . 

The Higgs self-coupling in the SM follows from equat- 
ing out the Higgs potential after the Higgs doublet is ex- 
panded around the electroweak symmetry breaking vac- 
uum expectation value, H = (0,t> + h) T '/y/2 in unitary 



gauge: 



2.2 



D -m h h 



(1) 



where 



"A related analysis of Higgs physics at a future linear collider can 
be found in Ref. El. 



r\v 2 )2, and v 2 = —^i 2 /rj. Since symmetry 
breaking in the SM relies on fi 2 < and r] > 0, the partial 
experimental reconstruction of the Higgs potential via 
the measurement of the trilinear Higgs vertex and its 
comparison to SM quantities (e.g. 2rj — g 2 m 2 ,/m 2 v ) is 
necessary to verify symmetry breaking due to a SM-like 
Higgs sector. Strictly speaking, a similar program needs 
to be carried out for the quartic Higgs vertex to fully 
reconstruct the Higgs potential by measurements. This 
task is, however, even more challenging due to an even 
smaller cross section of triple Higgs production [l~3|, [l4[ • 

A process at hadron colliders which is sensitive to 
the trilinear Higgs coupling is the previously mentioned 
Higgs pair production pp — > hh + X via gluon fusion. 
Consequently, this process has already been studied in 
the literature in detail [l4T - fr7| . From the known results 
it is clear that the LHC's potential towards measuring the 
trilinear coupling only in a single channel is insufficient. 
Combining different Higgs decay final states improves the 
situation, but exhausting the entire LHC search poten- 
tial requires also discussing the machine's capability to 
constrain the trilinear Higgs coupling in different kine- 
matical regimes. 

Different properties of signal and background pro- 
cesses in e.g. boosted final states as compared to un- 
boosted kinematics allows us to access Higgs decay chan- 
nels which are impossible to isolate from the background 
in a more inclusive search. This has impressively been 
demonstrated in the context of subjet-based analysis 
techniques [l8[ which have proven successful during the 
7 TeV LHC run In addition to that, initial state 

radiation can be an important effect when considering 
energetic final states. A dihiggs system recoiling against 
a hard hadronic jet accesses an entirely new kinemati- 
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FIG. 1: Sample Feynman graphs contributing to pp — > hh + X. Graphs of type (a) yield vanishing contributions due to color 
conservation. 



cal configuration^ , which is characterized by a large di- 
higgs invariant mass, but with a potentially smaller Higgs 
s-channel suppression than encountered in the back-to- 
back configuration of gg — > hh. 

The goal of this paper is to provide a comparative 
study of the prospects of the measurement of the trilinear 
Higgs coupling applying contemporary simulation and 
analysis techniques. In the light of recent LHC measure- 
ments, we focus our eventual analyses on mu = 125 GeV. 
However, we also put this particular mass into the con- 
text of a complete discussion of the sensitivity towards 
the trilinear Higgs coupling over the entire Higgs mass 
range rrih i5 1 TeV. As we will see, rrih — 125 GeV is a 
rather special case. Since Higgs self-coupling measure- 
ments involve end-of-lifctimc luminosities we base our 
analyses on a center-of-mass energy of 14 TeV. 

We begin with a discussion of some general aspects 
of double Higgs production, before we review inclusive 
searches for rrih — 125 GeV in the pp — > hh + X channel 
in Sec. Ill CI We discuss boosted Higgs final states inpp — > 
hh+X in Sec. Ill Dl before we discuss pp — > hh+j+X with 
the Higgses recoiling against a hard jet in Sec. Mil Doing 
so we investigate the potential sensitivity at the parton- 
and signal-level to define an analysis strategy before we 
apply it to the fully showered and hadronized final state. 
We give our conclusions in Sec. IIVI 



ggh and gghh interactions [2C 



4 3n p 

which upon expansion leads to 

LD+ 4 3^ G ^ G h 4 6W»" 



(2) 



(3) 



Studying these operators in the hh + X final state should 
in principle allow the Higgs self-coupling to be con- 
strained via the relative contribution of trilinear and 
quartic interactions to the integrated cross section. Note 
that the operators in Eq. (J3J have different signs which 
indicates important interference between the (nested) 
three- and four point contributions to pp — > hh + X al- 
ready at the effective theory level. 

On the other hand, it is known that the effective theory 
of Eq. ^ insufficiently reproduces all kinematical prop- 
erties of the full theory if the interactions are probed 



at momentum transfers Q 



2 > 



llll and the massive 



quark loops are resolved. Since our analysis partly re- 
lies on boosted final states, we need to take into account 
the full one-loop contribution to dihiggs production to 
realistically model the phenomenology. 



II. HIGGS PAIR PRODUCTION AT THE LHC 
A. General Remarks 

Inclusive Higgs p air production has already been stud- 
ied in Refs. [14T - [l7| so we limit ourselves to the details 
that are relevant for our analysis. 

Higgs pairs are produced at hadron colliders such as 
the LHC via a range of partonic subprocesses, the most 
dominant of which are depicted in Fig.Q] An approxima- 
tion which is often employed in phenomenological studies 
is the heavy top quark limit, which gives rise to effective 



B. Parton-level considerations 

In order to properly take into account the full dynam- 
ics of Higgs pair production in the SM we have imple- 
mented the matrix element that follows from Fig. [T] in 
the Vbfnlo framework [2l| with the help of the Fey- 
nArts/FormCalc/LoopTools packages [22|, with 
modifications such to include a non-SM trilinear Higgs 
coupling^. Our setup allows us to obtain event files ac- 
cording to the Les Houches standard [23], which can be 
straightforwardly interfaced to parton showers. Decay 
correlations are trivially incorporated due to the spin-0 
nature of the SM Higgs boson. 



^The phenomenology of such configurations can also be treated sep- 
arately from radiative correction contributions to pp — > hh + X. 



tThe signal Monte Carlo code underlying this study is planned to 
become part of the next update of Vbfnlo and is available upon 
request until then. 
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FIG. 2: Comparison of the (normalized) pt,h distributions in pp — > hh + X at LO for different multiples of the trilinear Higgs 
coupling A (m t = 172.5 GeV and mj = 4.5 GeV using CTEQ611 parton densities). 
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FIG. 3: Comparison of pp — > hh + X at LO. We choose 
m t — 175 GeV as in Ref. [l5|], from which we also obtain 
the dashed blue reference line, and mj, = 4.5 GeV and we use 
the CTEQ6U parton distributions. 

The resulting inclusive hadronic cross sections are plot- 
ted in Fig. [3j where we also show results for non-SM tri- 
linear couplings, varied around the SM value (see Eq. ([1])) 




Note that choosing a value different from Asm does not 
yield a meaningful potential in terms of Eq. ([T]), but al- 
lows to constrain A in hypothesis tests using, e.g., the 
CL S method [p. 

We also show the result of Ref. [l5| for comparison 
and find excellent agreement in total, keeping in mind 
that the results of Ref. [lj| were obtained using the GRV 
parametrizations of parton luminosities [25| . which are 
different from the CTEQ611 26;j set that we employ for 
the remainder of this paper§. 



Interference between the different non-zero contribu- 
tions depicted in Fig. [3J becomes obvious for the differ- 
ently chosen Higgs self-couplings. We also learn from 
Fig. [3J that the dihiggs cross section has a fairly large 
dependence on the particular value of the trilinear cou- 
pling for a to/j = 125 GeV Higgs boson. The qualitative 
Higgs mass dependence for different values of the trilinear 
self-coupling in Fig. [3J is easy to understand: The Higgs 
propagator in Fig.Q](c) is always probed off-shell at fairly 
large invariant masses; this renders the triangle contribu- 
tions subdominant compared to the box contributions of 
Fig. Q] (b) . For Higgs masses close to the mass of the loop- 
dominating top quark, we have s ~ 4mj, which results 
in resonant contributions of the three-point functions of 
Fig. HHc), well-known from one-loop gg — » h produc- 
tion [27|. This ameliorates the s-channel suppression of 
the trilinear coupling-sensitive triangle graphs and causes 
the dependence of the cross section on the trilinear cou- 
pling to become large at around < m*. 

To gain sensitivity beyond total event counts, it is im- 
portant to isolate the region of phase space which is most 
sensitive to modifications of the trilinear coupling in or- 
der to set up an analysis strategy which targets the tri- 
linear self-coupling most effectively. At the parton level, 
there is only a single phenomenologically relevant observ- 
able to hh production, which can be chosen as the Higgs 
transverse momentum pr,h- In Fig. [5] we show the dif- 
ferential pr,h distribution for different values of A and 
to/i = 125 GeV. The dip structure for A > Asm results 
again from phase space regions characterized by s ~ 4toj , 
which are available if < m t , and the resulting maxi- 
mally destructive interference with the box contributions. 

The above points suffice to give a qualitative assess- 
ment of the prospects of measurements of A in the pp — ¥ 
hh + X channel: 

• the Higgs bosons from inclusive dihiggs productions 



§ Using the integration-mode of FormCalc/LoopTools with the 



CTEQ611 set we obtain perfect agreement. 
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£=0£=1£=2 


bbWW 


ratio to £ 


= 1 


cross section before cuts 


59.48 28.34 13.36 


877500 


3.2 • 10" 


-5 


1 isolated lepton 


7.96 3.76 1.74 


254897 


1.5 • 10" 


-5 


MET + jet cuts 


1.54 0.85 0.44 


66595.7 


1.2 • 10" 


-5 


hadronic W reconstruction 


0.59 0.33 0.17 


38153.3 


0.9- 10" 


-5 


kinematic Higgs reconstruction 


0.028 0.017 0.007 


205.1 


8.3 • 10" 


-5 



TABLE I: Signal and background cross sections in fb for hh — > bbW + W . The Higgs self-coupling is scaled in multiples of the 
Standard Model value A = £ x Asm, Eq. @. The background is bbW + W~ production discussed at NLO in Ref. [2£| (K ~ 1.5). 



are naturally boosted pr,h ^ 100 GeV, 

• interference leads to an a priori A-sensitive phe- 
nomenology for rn /j ~ 125 GeV, 

• identical interference effects also cause the bulk of 
the sensitivity to A to follow from configurations 
with pT,h ~ 100 GeV, while the px,h shape at large 
values becomes similar for different values of A due 
to decoupling the triangle contributions at large 
partonic vl, 

• the cross section shows a dependence on the tri- 
linear coupling of Aa/asM — 50% when varying 
Ae [0,2A SM ]. 

We conclude our parton-level discussion of dihiggs pro- 
duction by noting that the higher-order corrections [16j . 
which result in a total K factor of ct nlo /ct lo > 1.85, re- 
sult from a large contribution from real parton emission. 
This is a characteristic trait of processes involving color- 
singlet final states at leading order, for which plenty of 
phase space for extra parton emission in addition to new 
initial-state parton combinatorics becomes available at 
next-to- leading order (NLO). Similar observations have 
been made for pp -t VV + X, where V = W ± ,Z, 7 [29j . 
The discussed characteristics of pp — > hh + X are there- 
fore not distorted when including NLO precision, and the 
parton shower will capture the characteristic features of 
the cross section upon normalizing to the NLO inclusive 
rate. 



C. Inclusive Higgs Pair Searches 

To measure A for a 125 GeV Higgs we need to iso- 
late modifications of pp — > hh + X cross sections around 
a NLO (hh + X) = 28.4 fb from the Higgses' exclusive de- 
cay channels. Given the small total inclusive cross sec- 
tion^, it is clear that even for s/s = 14 TeV and a tar- 
get luminosity of 0(lOOOfb _1 ) we need to focus on the 
Higgs decay channels with the largest branching ratios to 
phenomenologically visible final states to observe pp — > 



^The total inclusive single Higgs production cross section is 
16.5 pb [3dj for comparison. 



hh + X. These are [H h bb (59.48%), h -> W+W+ 
(20.78%), and h ->■ tt (6.12%). The decay h ^ ZZ 
(2.55%) is limited by the decays of the Z bosons to 
the clean leptonic final states Z — > e + e~,{i + ii~ (6.67%) 
(yielding BR(h -> e+e'+^+pr) = 0.013%). Hadronic Z 
decay modes can only be accessed in the boosted regime, 
which is not feasible for mh — 125 GeV [Hj]. 

We do not consider the final state hh — > 6677. A fea- 
sibility study for this particular channel was already pre- 
sented in Ref. @ . A realistic assessment of the sensitivity 
in 6677 depends on a realistic simulation of the diphoton 
fake rate due to multijet production, which is the dom- 
inant background to such an analysis, similar to Higgs 
searches in h — > 77. Details of the photon identification 
rely on the detector properties and the event selection 
approach, and we cannot address these issues in a realis- 
tic fashion. We focus in on h — > bb, W + W~ , t + t~ in the 
following. 

We generate all (fully showered and hadronized) back- 
ground samples with Sherpa [111 or MadEvent [34| . 
The signal events are interfaced to HerwigH — h [35| for 
showering and hadronization. 

1. hh-> bbbb 

The exclusive decay to two bb pairs is the most obvious 
channel to check for sensitivity due to its large branching 
ratio for rrih = 125 GeV. Using 6-tagging, it also possi- 
ble to access the intermediate invariant Higgs masses, 
for which the modifications due to A ^ Asm are well- 
pronounced. 

Passing the trigger-level cuts is not a problem for the 
signal events: the Higgses are naturally boosted and the 
PT-ordered b jets typically pass the staggered cuts on the 
transverse momentum pr,ji > 100 GeV, pr,j 2 > 80 GeV, 
Ptj 3 > 50 GeV, pr,j 4 > 40 GeV. However, as already 
mentioned, there is only one relevant scale to dihiggs pro- 
duction, and therefore our options to compete with the 
gigantic QCD pp —> bbbb + X background are highly lim- 
ited. Note that both Higgs bosons need to be recon- 
structed in order to be sensitive to the modifications of 
the trilinear coupling. 

In total, inclusive dihiggs production with decay to 
four b quarks has a signal-over-background ratio S/B 
which is too bad to be a suitable search channel, al- 
ready when we focus only on the QCD-induced four b 
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£ = £ = U = 2 


bbbb [QCD] bbbb [ELW] 


bbbb [QCD/ELW] 


ratio to £ 


= 1 


cross section before cuts 


59.48 28.42 13.36 


21165 


16.5 


160.35 


1.3- 10" 


3 


trigger+no leptons 


17.93 10.21 5.31 


5581.2 


8.0 


38.85 


1.8- 10" 


-3 


fat jet cuts 


13.73 8.23 4.50 


4761.0 


7.50 


31.65 


1.7- 10" 


-3 


first Higgs rec + 2b 


1.55 1.02 0.60 


235.22 


0.75 


1.32 


4.2 ■ 10" 


-3 


second Higgs rec + 2b 


0.137 0.094 0.059 


9.72 


0.011 


0.050 


9.6 • 10" 


-3 



TABLE II: Signal and background cross sections in fb for hh bbbb for boosted kinematics. The Higgs self-coupling is scaled 
in multiples of the Standard Model value A = £ x Asm, Eq. ((4]). Signal and background are normalized to the respective 
NLO cross sections. NLO bbbb production has been provided in Ref. [36| (inclusive K ~ 1.5). The mixed QCD+electroweak 
production and the purely electroweak production is currently not know at NLO QCD precision, and we therefore use an 
identical correction factor as for the QCD-induced production. 



background. Hence, this channel is not promising at the 
inclusive level and we revise it for boosted kinematics in 
Sec.QTDl 



too large for these cuts to bring down the S/B for this 
channel to a level for it to be useful in constraining the 
Higgs trilinear coupling. 



hh ->■ bbW+W- 



D. Boosted Higgs Pair Searches 



To gain multiple phenomenological handles to deal 
with the contributing backgrounds while preserving a sig- 
nal rate as large as possible we focus on h — » bb and 
h — » (W — » jj){W — > vl). The contributing background 
processes are tt and bbW + W~ production and we can 
employ cuts on missing transverse energy (MET) , lepton 
identification and px, and the reconstructed W resonance 
to reduce them. 

We require exactly one isolated lepton with pj>,i > 
10 GeV in the central part of the detector \y\ < 2.5, 
where isolation means an hadronic energy deposit of 
Er,had < 0.1 Et,i within a cone of R = 0.3 around the 
identified light-flavor lepton. In the next step we recon- 
struct the missing transverse energy (MET) $ T from 
all visible final state objects within the rapidity cover- 
age \y\ < 4.5 and require $ T > 20 GeV. Then, we 
use the anti-kT [37| algorithm as implemented in Fast- 
Jet [38[ (which we use throughout this paper) to re- 
construct jets with R = 0.6 and pr > 40 GeV, and 
require at least four jets in \y\ < 4.5. Afterwards we 
reconstruct the W boson by looping over all jet pairs. 
The jet pair that reconstructs the W mass best within 
60 GeV < rrijj < 100 GeV is identified as the W boson, 
and we subsequently remove these jets from the event. 
Analogous to the W reconstruction we reconstruct the 
Higgs within 110 GeV < mjj < 140 GeV. To reduce 
the backgrounds and identify the signal contributions we 
use a double 6-tag for the two jets which reconstruct the 
Higgs best. We use an efficiency of 60% with a fake rate 
of 2% in \y\ < 2.5 [39]. Thus, if one of the Higgs jets is 
outside \y\ < 2.5 our reconstruction fails. 

The results of this analysis flow can be found in Tab.Q] 
While the cuts bring down the background by a factor of 
4 x 10 3 , they also reduce the signal by nearly the same 
amount (1.5x 10 3 ). The requirement of two reconstructed 
Higgses has a strong effect on the background, however 
the initial cross-section (inclusively generated) is simply 



Moving on to the discussion of boosted final states, we 
can potentially gain sensitivity in the dominant Higgs de- 
cay modes, i.e. in the bb channels [l8[. The downside, 
however, is that we lose sensitivity to modifications of the 
trilinear couplings for harder Higgs bosons along the lines 
of Sec. Ill Bl Nonetheless, a measurement of the magni- 
tude of the dihiggs cross section is already an important 
task in itself. 

Recently, in the so-called BDRS analysis [l8j], it 
has been shown that applying jet substructure tech- 
niques on fatjets is a powerful tool to discriminate 
boosted electroweak-scale resonances from large QCD 
backgrounds. The BDRS approach proposes to re- 
combine jets using the Cambridge- Aachen (C/A) algo- 
rithm [42[ IHI with a large cone size to capture all decay 
products of the boosted resonance. Then one works back- 
ward through the jet clustering and stops when the clus- 
tering meets a so-called "mass-drop" condition: m,-, < 



31 

3 "31.32 

is not met 



VmiAR 



> 



lirrij witn /i = u.oO ana mm(j?y j l iPr j 2 
y CVL t using y cut = 0.09. If this condition 
the softer subjet fa is removed and the subjets of j\ 
are tested for a mass drop. As soon as this condition 
is met browsing backward through the cluster history 
the algorithm stops. In a step called "filtering" the con- 
stituents of the two subjets which meet the mass drop 
condition are recombined using the (C / A) algorithm with 
i?Bit = min(0.3, R b i/2). Only the three hardest filtered 
subjets are kept to reconstruct the Higgs boson and the 
two hardest filtered subjets are ^-tagged. The filtering 
step reduces the active area of the jet tremendously and 
makes the Higgs-mass reconstruction largely insensitive 
to underlying event and pileup. 

For the reconstruction of the boosted Higgs bosons 
in Sec. HID II and Sec. HID 21 we adopt this approach 
without modifications. It is worth noting that other 
techniques, possibly in combination with the BDRS ap- 
proach, can improve on the Higgs reconstruction em- 
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c = o e = i £ = 2 


bbrr bbrr [ELW] 


bbW + W~ 


ratio to £ 


= 1 


cross section before cuts 


59.48 28.34 13.36 


67.48 


8.73 


873000 


3.2 ■ 10" 


-5 


reconstructed Higgs from rs 


4.05 1.94 0.91 


2.51 


1.10 


1507.99 


1.9 ■ 10" 


-3 


fatjet cuts 


2.27 1.09 0.65 


1.29 


0.84 


223.21 


4.8 • 10" 


-3 


kinematic Higgs reconstruction (m b i) 


0.41 0.26 0.15 


0.104 


0.047 


9.50 


2.3 • 10" 


-2 


Higgs with double 6-tag 


0.148 0.095 0.053 


0.028 


0.020 


0.15 


0.48 





TABLE III: Signal and background cross sections in fb for hh — > bbr + T~ for boosted kinematics. The Higgs self-coupling is 
scaled in multiples of the Standard Model value A = £ x Asm, Eq. ((4]). The background comprises ti with decays to t — ¥ brv T , 
and bbr + T^ for pure electroweak and mixed QCD-electroweak production, normalized to the respective NLO rates. The 
bbW + W~ NLO cross sections are provided in [23] (K ~ 1.5), for the mixed and the purely electroweak contributions we infer 
the corrections from Zbb (K ~ 1.4) and ZZ (K ~ 1.6) production using mcfm msn. 



ciency and can therefore increase the sensitivity of the 
following searches 44 1. 



1. hh-^ bbbb 



As already pointed out, the Higgs bosons are natu- 
rally boosted, and requiring two fatjets subject to BDRS 
tagging [l8T ] can improve the very bad S/B in the con- 
ventional pp — > bbbb + X search without losing too much 
of the dihiggs signal cross section. 

In the analysis, we veto events with light leptons 
Pt,i > 10 GeV in \y\ < 2.5 to reduce ti, where the 
leptons are again assumed isolated if £r,had < 0.1Et,i 
within R < 0.3. We need to make sure that the events 
we want to isolate pass the trigger level. For this reason, 
we recombine final state hadrons to jets with R = 0.4 and 
Pt > 40 GeV and require at least four jets and the fol- 
lowing staggered cuts: pT,j x > 100 GeV, pr,j 2 > 70 GeV, 
Ptj 3 > 50 GeV. All jets have to be within detector cov- 
erage \y\ < 4.5. 

For the events that pass the trigger cuts, we apply 
a "fatjet" analysis, i.e. require at least two jets with 
Ptj > 150 GeV and R — 1.5 in the event. We apply 
the BDRS approach to both of these fatjets using /i cut = 
0.66 and y cut = 0.09. The reconstructed Higgs jets need 
to reproduce the Higgs mass within a 20 GeV window: 
115 GeV < mh < 135 GeV, and we additionally require 
that the two hardest filtered subjets are 6-tagged. 

We generate the backgrounds with exclusive cuts to 
make our cut-analysis efficient, yet inclusive enough to 
avoid a bias. More precisely we demand two pairs of b 
quarks to obey R bb < 1.5, pr(bb) > 100 GeV, m(bb) > 
50 GeV, while p T . b > 20 GeV, and \y b \ < 2.5. The 
(anti-)fc are generated with R bb > 0.2. 

The results are collected in Tab. |TTJ Again, while the 
cuts allow an improvement in S/ B by an nearly an order 
of magnitude, we are still left with a small signal rate on 
top of a very large background so that this channel is in 
the end also not promising. 



2. hh -> bbr+r- 

A promising channel is dihiggs production with one 
Higgs decaying to a pair of r leptons. This decay chan- 
nel in association with two jets is one of the main search 
channels for single light Higgs production [53,[48[ and has 
recently been used to put bounds on Higgs production 
by Cms [49j]. The reconstruction of r leptons is delicate 
from an experimental point of view, and current analysis 
strategies mostly rely on semi-hadronic r pair decays in 
the context of Higgs searches (see e.g. Ref. [4j|). The r 
identification is performed using likelihood methods 
which do not allow a straightforward interpretation in 
terms of rectangular cuts used in e.g. Ref. 48]. Con- 
sequently, with likelihood r taggers unavailable to the 
public, a reliable and realistic estimate is hard to obtain. 
For this reason, we choose a r reconstruction efficiency of 
80% with a negligible fake rate. This is not too optimistic 
in the light of the likelihood approaches of Ref. [50| , bear- 
ing in mind that our analyses are based on end-of-lifetime 
luminosities, for which we may expect a significant im- 
proved r reconstruction when data is better understood. 
We choose a large enough Higgs mass window for the 
reconstruction, in order to avoid a too larg e systematic 
pollution due to our assumption (in Ref. [49[ CMS quotes 
a O(20%) of the reconstructed Higgs mass). 

In more detail, we require two r jets with pt > 20 GeV, 
reproducing the Higgs mass within 50 GeV, m TT = nth ± 
25 GeV. Then we use the C/A algorithm to reconstruct 
fatjets with R — 1.5 and pxj > 150 GeV and require 
at least one fatjet in the event. Thereby we demand 
the fatjets to be sufficiently isolated from the rs. We 
subsequently apply the BDRS approach to the fatjet with 
Mcut = 0.66 and y cut = 0.09. The two hardest filtered 
subjets need to pass b tags and the reconstructed Higgs 
jet has to be in mh ± 10 GeV. £>-tagging is performed for 
\y\ < 2.5 and we assume an efficiency of 70% and a fake 
rate of 1% following Ref. 51| . 

We generate the bbrr and pure electroweak bbrr back- 
grounds with exclusive cuts to make our cut-analysis rea- 
sonably efficient , yet inclusive enough to avoid a bias. 
More precisely we demand the two b quarks to obey 
R bb < 1.5, p T (bb : TT) > 100 GeV, to(66,tt) > 50 GeV, 
while pr,b.T > 20 GeV, and \y b . T \ < 2.5. The bs and 
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FIG. 4: Sample Feynman graphs contributing to pp — > hh + j + X. Not shown are the qg, qg and qq subprocesses. 



b 



effective theory 
full theory 



0.1 



0.01 r 



0.001 




200 300 
Ptj [GeV] 



r,oi) 



10 




200 300 
p T ,j [GeV] 



500 



FIG. 5: Comparison of the leading order prj spectrum for 
pp —¥ hh + j + X production. Shown are distributions for the 
effective interaction (obtained with MadGraph v5 [34| via 
FeynRules and Ufo HI]), and the full one- loop matrix 
element calculation. We again choose rat = 172.5 GeV and 
nib = 4.5 GeV using CTEQ611 parton densities and /if = 
fJ-R = PT,j + 2m h . 



rs are generated with Rbb,rr > 0.2. On the other hand, 
the bbW~ W + sample is generated inclusively, and is the 
same sample used in the unboosted bbW~W + analysis 
in the previous section. 

The results are shown in Tab. IIIII The initial back- 
ground cross-section looks very large due to it being in- 
clusively generated. However, once we take into account 
the small branching ratio of W — > tv this drops dramati- 
cally. After requiring two 6-tagged jets which reconstruct 
the Higgs mass we are left with an S/B of nearly half for 
the £ = 1 case (and nearly one in for £ = 0). The cross- 
section is also reasonable, corresponding to 95 events for 
1000 inverse femtobarns of luminosity. This channel is 
hence very promising indeed. 



III. HIGGS PAIR PRODUCTION IN 
ASSOCIATION WITH A HARD HADRONIC JET 

A. Parton-Level considerations 

The qualitatively poor agreement of the effective the- 
ory of Eq. ([3]) with the full theory persists if additional 
jet radiation is included. Naively we could have expected 
that accessing smaller invariant masses in the Higgs sys- 
tem due to significant initial state radiation might re- 
sult in a better agreement with the effective theory of 
Eq. ([3]). However, especially for hard jet emission, which 
allows the Higgs pairs to access large invariant masses 
in a new collinear kinematical configuration compared 
to pp — > hh + X, the disagreement of full and effective 
theories is large (Fig. [SJ. 

Given these shortcomings of the effective theory, we 
implement the full matrix element in the Vbfnlo frame- 
work using FeynArts/FormCalc/LoopTools. We 
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FIG. 6: Comparison of the (normalized) leading order maxp^j, distributions in pp — s- hh + j + X for different multiples of the 
trilinear Higgs coupling A (m« = 172.5 GeV and mj = 4.5 GeV using CTEQ611 parton densities), and pr,j > 20 (100) GeV in 
the upper (lower) row, respectively. Factorization and renormalization scales are chosen /if = fJ-R = Pr.j + irrih- 
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FIG. 8: Comparison of leading order pp — > hh + j + X production cross sections for three different values of the Higgs 
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have checked our phase space implementation for the 
effective theory's matrix element against MadEvent. 
Some of the contributing Feynman graphs to the dom- 
inant gg-initiated subprocess are shown in Fig. @] note 
that again only a subset of the contributing diagrams 
is sensitive to non-standard hhh couplings. Interference 
between these and the remaining contributions is again 
obvious from Fig. [5] especially at around < m t , which 
can again be explained along the lines of Sec. Ill B[ 

In comparison to pp — > hh + X , we find sizably larger 
dependence on A of the total cross section, Fig [5] For 
Ptj > 20 GeV we have Aa/a SM ~ 100% for a variation 
< A < 2Asm- This is due to the larger available phase 
space for the dihiggs system. The intermediate s chan- 
nel Higgs in Fig. 0] (a), (c) is probed at smaller values 
compared to Fig. [T] (c), suppression is ameliorated and 
(destructive) interference becomes more pronounced. 

With a dihiggs system that becomes less back-to- 
back for increasingly harder jet emission, the character- 
istic dip structure encountered in the px,h spectrum of 
pp —> hh + X is washed out (Fig. [5]). Characteristic 
imprints can still be observed in the dihiggs invariant 
mass or, equivalently, in the dihiggs separation in the 
azimuthal-angle — pseudorapidity plane, Fig. [7] 

Let us summarize a few points relevant to the analysis 
of pp — > bbbb + j + X production for ptj > 100 GeV: 

• The dihiggs+jet cross section has a comparably 
large dependence on the value of the trilinear cou- 
plings as compared to pp — > hh + X (Act/ctsm — 
45% when varying A € [0, 2Asm]), 

• the sensitivity to non-standard values of the trilin- 
ear coupling arises from phase space configurations 
where the two Higgs bosons are close to each other 
in the central part of the detector, i.e. for rather 
small values of the invariant masses, 



• as a consequence, the hadronic Higgs decay prod- 
ucts are likely to overlap, and to fully reconstruct 
the busy hh decay system we need to rely on jet- 
substructure techniques. 



Let us again comment on the impact of higher order 
QCD contributions. A full NLO QCD computation for 
pp — > hh+j+X is yet missing, but most pp —¥ VV+j+X 
(V = W ± ,Z,'-f) production cross sections, which have 
similar properties from a QCD point of view, are known 
to NLO QCD precision [52]. Also, the NLO QCD cross 
sections ior pp — > Vh+j+X (V = W ± 1 Z) have been pro- 
vided in Ref. [Hj]. Given that the QCD sector is largely 
agnostic about the matrix elements' precise electroweak 
properties (taken apart the partonic composition of the 
initial state) , it is not a big surprise that all these produc- 
tion cross sections exhibit a rather similar phenomenol- 
ogy at NLO QCD. The total inclusive K factors range 
around K ~ 1.3 and result from unsuppressed parton 
emission. It is hence reasonable to expect the QCD cor- 
rections to pp — > hh + j + X to be of similar size, and 
parton shower Monte Carlo programs to reasonably re- 
produce the dominant kinematical properties. 

For the remainder of this paper we do not include the 
weak boson fusion component [fj4| to one-jet inclusive 
production. This the second largest contribution to inclu- 
sive dihiggs production, but it is still smaller than hh + j 
production from gluon fusion in the phase space region 
we are interested in. For maxprj > 80 GeV we have 
&WBF(hh + 2j) ~ 0.5 fb in the SM, so this amounts to a 
O(+10%) correction to our inclusive signal estimate (well 
inside the perturbative uncertainty of pp — > hh + j + X). 
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e = (H = ie = 2 


bbbbj [QCD] bbbbj 


[QCD/ELW] bbbbj [ELW] 


ratio to £ 


= 1 


cross section before cuts 




6.45 3.24 1.81 


29400 


513.36 


10.0 


1.1 • 10" 


4 


trigger+fatjet cuts 




1.82 1.08 0.69 


10579.8 


211.04 


4.16 


1.0 • 10" 


-4 


first kinematic Higgs rec (new tag 


ger) + 2b 


0.30 0.20 0.13 


331.84 


10.82 


0.54 


0.5 • 10" 


3 


sec kinematic Higgs rec + 2b 




0.09 0.059 0.039 


54.1 


2.46 


0.066 


1.0 ■ 10" 


3 


invariant mass + pr,j cut 




0.049 0.031 0.022 


36.06 


0.92 


0.030 


0.9 • 10" 


3 



TABLE IV: Signal and background cross sections in fb for hh+j^ bbbb + j for boosted kinematics. The Higgs self-coupling 
is scaled in multiples of the Standard Model value A = £ x Asm, Eq. ((4]). None of the contributing backgrounds' normalization 
is known to NLO QCD precision. We therefore include a conservative correction factor of K — 2. 



B. Boosted Higgs searches in association with a jet 

From Fig. [6] we see that the Higgs bosons are again 
naturally boosted. Events of this signature possess an 
hadronically more active final state. The pure BDRS ap- 
proach works very well if there is no other hard radiation 
inside the fatjet except the decay products of the Higgs 
boson. Here it is likely that the additional hard jet ends 
up in one of the fatjets challenging a good reconstruction 
of the Higgs. Therefore, we modify the tagger similar 
to the Higgs tagger in [55] : The last clustering of the 
jet j is undone, giving two subjets ordered such 

that m,j 1 > nij 2 . If nij 1 > 0.8rrij we discard j% and keep 
ji, otherwise both j\ and ji are kept. For each subjet 
ji that is kept, we either add it to the list of relevant 
substructures (if m, < 30 GeV) or further decompose it 
recursively. After performing this dcclustering procedure 
(we do not stop after observing a mass drop but continue 
to decluster the jets until we obtain a set of hard sub- 
jets inside the fat jet) we recombine the constituents of 
every two-subjet combination with the C/A algorithm 
using i?fiu = min(0.3, Ai^jj j2 /2). For each combination 
we keep the three hardest filtered subjets and call it a 
Higgs candidate. The two hardest filtered subjets of the 
Higgs candidate with the mass closest to the true Higgs 
mass of 125 GeV we require to be 6-tagged. We find 
that this tagger recovers roughly 40% more of the signal 
events while keeping S/B constant. 



1. hh + j ->• bbbb + j 

We proceed in a similar manner to the analysis outlined 
in Sec. Ill Dl but with modifications of the Higgs tagger 
in order to preserve a larger signal cross section. 

We generate the backgrounds with the following 
parton-level cuts, yet inclusive enough to avoid a bias. 
We require that two bb combinations obey pr{bb) > 
100 GeV and m(bb) > 100 GeV, while |y 6 | < 2.5 and 
PT,b > 20 GeV. The 6s are separated by R bb > 0.2. The 
additional jet is generated with pr > 80 GeV in \yj\ < 4.5 
and is separated from the bs by AR > 0.7. Signal events 
are generated with ptj > 80 GeV. 

In the analysis, we again veto events with isolated lep- 
tons in \y\ < 2.5 for p T j > 10 GeV and £r,had < 0.1 E T j 



within R < 0.3. To assure the trigger requirements 
we ask for at least five jets and the following staggered 
cuts on the transverse momentum: prji > 120 GeV, 
Pt, J2 > 100 GeV, p T ,j 3 > 70 GeV, p T ,j 4 > 40 GeV. 
All jets have to fall inside the detector \y\ < 4.5. The 
events which pass these trigger criteria are again ana- 
lyzed in a subjet approach: We reconstruct fatjets with 
Pr.j > 150 GeV and R = 1.5. At least one fatjet 
has to be present which fulfills the mass-drop condi- 
tion [18j with an invariant mass of rrij > 110 GeV. The 
reconstructed Higgs has to be double 6-tagged and have 
p T ,h > 150 GeV. 

The hardest of these fatjets is declustered with a tag- 
ger which is inspired by the Higgs tagger of Ref. [55j . 
The reconstructed Higgs has to be double 6-tagged with 
Ptji > 150 GeV. We subsequently remove the recon- 
structed Higgs' constituents from the event. 

The remaining constituents are reclustered using the 
anti-kT algorithm with R = 0.4 and p T ,j > 30 GeV. The 
two jets which reconstruct the Higgs mass best within 
nih = 125 ± 10 GeV are ^-tagged with 60% tagging effi- 
ciency and 2% fake rate; 6-tagging is again performed in 
\y\ < 2.5. Then, the two jets are again removed from the 
event. For the two reconstructed Higgs jets we require 
an invariant mass (j>h 1 + Ph 2 ) 2 > 400 2 GeV 2 and the 
last remaining jet needs to be hard pr > 100 GeV. In 
total, this corresponds to a signal signature discussed in 
Sec. HO] 

The results of this analysis flow can be found in 
Tab. IIVI For the numbers quoted there it is clear that 
our cut setup serves to highly reduce the contributing 
backgrounds, but the QCD-induced cross sections have 
a too large initial value. In the QCD background, while 
the bbbb was predominantly gluon initiated, the bbbb + j 
receives large contributions from qg initial states, leading 
to final states with a large invariant mass. This in turn in- 
creases the amount of background in searches for boosted 
resonances. The QCD backgrounds can be reduced by a 
factor 0(1000) while the signal rate is decreased by a 
factor ~ 100. In total, this does not leave a large enough 
S/B to be relevant from the point of systematics, and we 
conclude that pp — > hh + j — > bbbb + j is not a sensitive 
channel. 
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bbr + T j bbr 


-T~j [ELW] 


ttj 


ratio to £ = 1 


cross section before cuts 


6.45 3.24 1.81 


66.0 


1.67 


106.7 


1.9 ■ 10" 2 


2 rs 


0.44 0.22 0.12 


37.0 


0.94 


7.44 


4.8 ■ 10" 3 


Higgs rec. from taus + fatjet cuts 


0.29 0.16 0.10 


2.00 


0.150 


0.947 


5.1 • 10~ 2 


kinematic Higgs rec. 


0.07 0.04 0.02 


0.042 


0.018 


0.093 


0.26 


2b + hh invariant mass + pr,j cut 


0.010 0.006 0.004 


<0.0001 


0.0022 


0.0014 


1.54 



TABLE V: Signal and background cross sections in fb for hh + j — > bbr + T~ +j for boosted kinematics. The Higgs self-coupling 
is scaled in multiples of the Standard Model value A = £ x Asm, Eq. (Q. The QCD corrections to tt + j have been discussed 
in Ref. [5r3 ] (K ~ 1.1). For the pure electroweak production we take the results of as a reference value (K ~ 1.3). The 
corrections to mixed production are unknown and we conservatively use a total inclusive QCD correction K = 2. 



2. hhj ->■ bbr* 



t J 



IV. SUMMARY 



We follow closely the steps described in Sec. Ill D 2 1 and 
Sec. HUB H 

We generate the backgrounds with the following 
parton-level cuts to have a reasonably efficient analy- 
sis, yet inclusive enough to avoid a bias. We require 
p T (bb, tt) ^ 100 GeV and m(bb, tt) > 90 GeV (100 GeV 
in case of tt + j), while \yb, T \ < 2.5 and pr,b,r > 20 GeV. 
The 6s and rs are separated by Rbb.TT > 0.2. The addi- 
tional jet is generated with px > 80 GeV in \yj\ < 4.5 
and is separated from the bs by AR > 0.7. Signal events 
are generated with ptj > 80 GeV. 

We require exactly two r jets in an event in \y T \ < 2.5 
with pt > 20 GeV and assume an identification efficiency 
of 80% each. The rs have to reconstruct to an invariant 
mass of m/!±25 GeV. Then we use the C/A algorithm to 
reconstruct fatjets with R = 1.5 and px,j > 150 GeV and 
require at least 1 fatjet in the event which is sufficiently 
isolated from the rs. Then we apply the Higgs tagger 
described in Sec. HID I and require the reconstructed Higgs 
jet have a mass of ± 10 GeV and pt.h > 150 GeV. To 
suppress the large tt background we reject events where 
the invariant mass of the two reconstructed Higgs bosons 
is below 400 GeV. After removing the constituents of 
the reconstructed Higgs bosons from the final state we 
cluster the remaining final state constituents using the 
anti-kT algorithm R = 0.4 and pr,j > 30 GeV. Finally, 
we require at least one jet with pr > 100 GeV. 

We find that these cuts can suppress the backgrounds 
confidently as long as the r fake rate is sufficiently small. 
Due the large invariant mass of the final state, several 
high-py jets and possibly leptons from the r decays we 
expect that these events can be triggered on easily. The 
full analysis flow can be found in Tab. [V] The initial 
background contributions are significantly lower, as this 
final state does not have a dominant purely QCD-induced 
component. In total we end up with an estimate on 
S/B ~ 1.5. This means that with a target luminosity 
of 1000 fb -1 , constraints can be put on A in this channel. 



We have studied the prospects to constrain the trilinear 
Higgs coupling by direct measurements at the LHC in 
several channels, focussing on mu = 125 GeV. This is 
also the mass region which is preferred by electroweak 
precision data, and where we currently observe excesses 
in data both at the LHC and the Tevatron. Depending on 
the particular decay channel, we find a promising signal- 
to-background ratio at the price of a very small event 
rate. 

Higgs self-coupling measurements for a SM Higgs in 
this particular mass range are typically afflicted with 
large backgrounds, so that achieving maximal sensitivity 
requires the combination of as many channels as possi- 
ble. For dedicated selection cuts we obtain signal cross 
sections in Higgs pair production of the order of 0.01 to 
0.1 fb and measurements will therefore involve large data 
sets of the 14 TeV run with a good understanding of the 
involved experimental systematics. 

Searches for unboosted kinematics of the Higgs bosons 
do not allow any constraint on the trilinear coupling or 
total cross-section to be made. However, requiring the 
two Higgses to be boosted and applying subjet methods 
to boosted pp — > hh+X and pp — > hh+j+X production, 
we find a sensitive S/ B particularly for final states involv- 
ing decays into rs. A necessary condition for sensitivity 
in these channels is a sufficiently good r reconstruction, 
but more importantly, a small fake rate. Unfortunately, 
while boosting the Higgses increases S/B, it leads us into 
a region of phase space which lacks sensitivity to the tri- 
linear coupling. 

In addition to inclusive dihiggs production we find that 
dihiggs production in association with a hard jet shows 
an improved sensitivity to the trilinear Higgs coupling. 
However to exploit this scenario still requires the use of 
boosted techniques which require thorough evaluation on 
data. 

Assuming the efficiency for r-tagging and the hadronic 
Higgs reconstruction as outlined in this work are con- 
firmed using data, the 66r + r~ and 66r + r~ + j channels 
can be used to constrain the Higgs self-coupling in the 
SM at the LHC with a data set of several hundred in- 
verse femtobarns. The analysis strategies developed in 
this paper will also help to improve bounds on dihiggs 
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production in scenarios with strong electroweak symme- 
try breaking and related models, which also predict en- 
hanced dihiggs production cross sections [1, H3| . 
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